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FOREWORD 


The  computer-aided  preliminary  design  system  for  light 
antitank  weapons  described  in  this  report  was  developed  by 
Athena  Engineering  Company  under  contract  DAAH01-73-C-0654  with 
the  U.S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama.  The 
development  was  sponsored  and  technically  monitored  by  the 
Ground  Equipment  and  Materials  Directorate,  US  Army  Missile 
Research,  Development,  and  Engineering  Laboratory.  The  computer- 
aided  design  system^  referred  to  as  CADLAW,  is  operational  on 
the  MICOM  CDC  6600  computer  and  operates  interactively  through 
a Tektronix  4015  communication  terminal.  The  system  is  also 
operational  on  a UNIVAC  1108  system  used  by  Athena  Engineering 
Company.  < 

The  author  wishes  to  acknowledge  the  valuable  technical 
contributions  of  Athena  Engineering  Company  personnel  Mr.  Jeunes 
Dagen,  Dr.  James  L.  Hill,  and  Dr.  Howard  B.  Wilson  and  Mr. 

Richard  Eppes,  Ground  Equipment  and  Materials  Directorate. 

Mr.  Dagen  provided  comprehensive  general  programming  support 
including  developing  the  interactive  graphics  capability.  Dr. 

Hill  developed  the  CADLAW  component  which  calculates  the  missile 

i 

trajectory.  Dr.  Wilson  developed  the  CADLAW  recoilless  launcher 
component  aifd  provided  the  basic  optimization  program,  and  Mr. 


Eppes  provided  overall  technical  guidance  and  assisted  in 
making  the  CADLAW  program  operational  on  the  CDC  system. 

Since  the  CADLAW  program  is  interactive  it  is  largely 
system  dependent  due  to  the  overlay  structure.  For  this  rea- 
son no  listing  of  the  program  is  included  in  this  report. 

The  program  can  be  made  operational  on  other  systems  with 


minor  modifications. 
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1 . INTRODUCTION 


Conventional  preliminary  design  of  missile  systems 
requires  the  interaction  of  people  representing  various 
technical  and  managerial  disciplines  in  a process  which  is 
sometimes  vaguely  defined  and  which  is  almost  always  itera- 
tive in  nature.  The  advantages  of  incorporating  the  computer 
into  this  design  process  appear  numerous:  the  process  could 

be  .conducted  more  rapidly  with  less  human  effort,  the  inter- 
action of  various  diSitclplines  could  be  examined  more  easily 
and  precisely,  a greater  number  of  iterations  could  be  per- 
formed in  a given  period  of  time,  a more  efficient  "optimum 

I 

design"  could  be  selected,  etc.  The  developments  reported 
here  are  directed  toward  the  goal  of  utilizing  the  computer 
as  an  aid  in  the  preliminary  design  of  light  antitank 
weapons . 

The  conventional  preliminary  design  process  typically 
begins  by  having  the  designer  select,  on  the  basis  of  pre- 
vious experience,  a design  which  appears  to  meet  the  design 
specifications.  After  the  initial  design  has  been  selected 
it  is  passed  to  various  people  representing  a broad  spectrum 
of  technical  and  managerial  disciplines  where  the  pekrformance 
and  financial  costs  are  evaluated.  Detailed  evaluation  of 


the  initial  design  generally  reveals  the  need  for  several 
design  modifications  whose  impact  must  be  evaluated  through 
trade-off  studies.  Since  the  goals  of  high  strength  and  per- 
formance conflict  with  the  concurrent  goals  of  low  weight 
and  cost,  inevitably  compromises  must  be  made. 

The  performance  evaluations  and  trade-off  studies 
conducted  in  the  course  of  the  design  development  are  usually 
conducted  separately  and  involve  a significant  amount  of 
human  effort.  Also,  the  final  selection  of  a particular 
design  configuration  is  generally  not  obvious  and  involves 
the  judgment  of  the  designer  to  a large  extent.  Use  of  the 
computer  to  aid  in  this  design  process  appears  an  attractive 
goal. 

Although  computers  are  employed  to  a large  ektent  in 
the  performance  evaluation  and  trade-off  studies  in  the 
conventional  design  process,  they  are  not  typically  used  in 
the  decision-making  process  nor  in  the  interactive  demand 
mode.  The  CAD-E  program  sponsored  by  NASA  [1]  and  the  CAMS 
program  sponsored  by  the  Air  Force  [2]  represent  attempts 
to  utilize  the  computer  to  a greater  extent  in  the  decision- 
making process  of  design. 

{■ 

To  explore  the  feasibility  of  automating  the  pre- 
liminary design  process  for  a realistic  missile  system,  a 

I. 

light,  shoulder-fired,  antitank  weapon  system  was  selected 
for  application.  The  objective  of  the  developmental  effort 
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was  to  demonst^rate  the  feasibility  of  the  computer-aided 
preliminary  design  process  on  a realistic  weapon  system. 

No  attempt  was  made  to  be  unduly  precise  in  system  modeling. 
Instead,  reasonable  mathemati  ^al  models  of  all  systems  were 
developed  with  the  awareness  that  more  sophisticated  models 
could  be  substituted  for  developmental  models,  if  desired. 

The  basic  operation  of  the  system  would  still  be  valid. 

Two  versions  of  the  light  antitank  weapon  were 
selected  for  application,  a conventional  in-tube  burning 
rocket  system  and  a closed  combustion  chamber  recoilless 

I 

system.  Two  versiors  were  studied  as  a result  of  contenqpor- 
ary  Interest  by  the  U.  S.  Army  Missile  Command.  The  in-tube 
burning  rocket  system  has  been  the  subject  of  several  recent 

developmental  studies  [3-6].  The  closed  recoilless i system 

>1 

has  also  attracted  attention  due  to  its  apparent  lack  of 
signature  and  blast  effects  [7-11]  . The  object  of  the  pre- 
sent effort  was  not  to  develop  optimum  designs  for  these 
type  weapons  but  rather  to  develop  a computer-aided  design 
system  which  would  allow  such  optimum  design  studies  to  be 
performed.  i 

A significant  feature  of  the  CADLAW  system  is  the 
ability  to  select  an  "optimum  design"  automatically.  It  is 
not  clear  that  conventional  preliminary  design  processes 

I 

result  in  such  an  optimum  design.  Thus,  the  ability  to 
determine  such  an  optimum  design  automatically  is  significant. 
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The  feature  which  allows  a determination  of  the  optimum 
design  is  based  upon  an  optimization  method  known  as  the 
"flexible  tolerance"  algorithm  [12].  This  method  con- 
stitutes a particularly  versctile  optimization  method  which 
is  a member  of  a wide  class  of  such  methods  [13] . The 
method  involves  an  iterative  pattern  search  of  a multi- 
dimensional hypersurface  specified  by  the  designer  in  the 
form  of  an  objective  function.  The  flexible  tolerance 
algorithm  searches  the  hypersurface  defined  by  the  objective 
function  until  an  optimum,  maximum  or  minimum,  design 
point  is  reached.  Thus,  an  optimum  design  depends  on  the 
character  of  the  objective  function  and  constraints  on  the 
design  variables. 

The  in-tube  burning  rocket  system  and  recoilless 

I 

system  have  propulsion,  structures  and  weights,  aerodynamics, 
guidance  and  control,  and  trajectory  simulation  components. 
From  an  optimization  standpoint  the  propulsion  and  the  struc- 
ture and  weight  components  are  separate  from  the  aerodynamics, 
the  guidance  and  control,  and  the  trajectory  simulation  com- 
ponents. The  design  variables  associated  with  the  propulsion 
and  the  structures  and  weights  systems  are  selected  in  an  op- 
timum  manner  independent  of  the  aerodynamics,  the  guidance  and 
control,  and  the  trajectory  system  variables.  This  separation 
is  permitted  in  this  application  since  all  propulsion  is  accom- 
plished within  the  launch  tube  before  the  aerodynamic  and 
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guidance  characteristics  affect  performance.  Hence,'  one 
objective  function  is  needed  for  the  propulsion  phase  of  the 
system  performance  and  another  for  the  flight  phase. 

In  summary,  computer-aided  design  techniques  offer 
advantages  Impossible  to  obtain  with  conventional  preliminary 
design  and  analysis  techniques.  Computer  modeling  of  all 
system  components  and  their  interactions  permit  design 
iterations  to  be  evaluated  very  rapidly.  The  speed  of 
evaluation  permits  optimum  designs  to  be  selected  in  an  auto- 
m2i.ted  fashion.  Perhaps  the  greatest  advantage  of  the  auto- 
mated design  cycle  '.*.s  the  freedor  from  computational  drudgery 
afforded  the  designer  who  is  then  able  to  devote  more  of  his 
effort  to  creative  thinking.  Details  of  the  CADLAW  system 
are  contained  in  the  following  sections  of  this  report. 
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2.  GENERAL  DESIGN  SYSTEM  DESCRIPTION 


As  indicated  previously,  the  CADLAW  system  operates 
in  a "demand  mode"  and  requires  interaction  of  the  designer 
with  the  computer.  For  operation  on  the  MICOM  CDC-6600 
computer  the  computer  program  was  designed  to  operate  in 
separate  segments  so  that  the  OVERLAY  feature  of  CDC  com- 
puters could  be  employed.  Table  2.1  contains  a list  of  the 
various  overlays  in  the  CADLAW  system  and  gives  an  indication 
of  the  basic  program  structure. 


TABLE  2.1 

CADLAW  SYSTEM  COMPONENTS 


OVERLAY 

(0,0) 

- 

Main  Program 

OVERLAY 

(1,0) 

- 

Conventional  Rocket  System 

OVERLAY 

(1,1) 

- 

Optimization  Algorithm 

OVERLAY 

(1,2) 

- 

Performance  Evaluation 

OVERLAY 

(2,0) 

- 

Recoilless  Rocket  System 

OVERLAY 

(2,1) 

- 

Optimization  Algorithm 

OVERLAY 

(2,?) 

- 

Performance  Evaluation 

OVERLAY 

(3,)) 

Trajectory  Simulation 

OVERLAY 

(3,1) 

- 

Optimization  Algorithm 

OVERLAY 

(3,-2) 

- 

Performance  Evaluation 

8 


4'i 


m- 


2.1  System  -Components 

The  Main  Program,  OVERLAY  (0,0)  serves  mainly  as  an 

Instructional  element  to  explain  the  basic  purpose  and  use 

i 

of  the  design  system.  The  user  is  provided,  interactively, 
with  information  about  the  system  and  is  allowed  to  select 
which  system  component  he  desires  to  work  with. 

It  is  assumed  that  all  propulsion,  in  both  the  con- 
ventional and  recoilless  systems, is  accomplished  within  the 
launch  tube.  This  assumption  allows  the  trajectory  Simula- 

I 

tion  and  guidance  and  control  phases  of  the  missile  perform- 
ance to  be  separateu  from  the  propulsion  phase.  Since  the 
strength  and  weight  requirements  are  largely  associated  with 
the  forces  exerted  during  the  propulsion  phase,  these  com- 
ponents are  also  separated  from  the  trajectory  simulation 
and  guidance  and  control  phases  of  the  system  performance. 

The  optimized  design  and  performance  evaluation  of 
the  conventional  in-tube  burning  light  antitank  propulsion 
and  launch  systems  is  accomplished  in  OVERLAY  (1,0)  of 
CADLAW.  The  influence  of  such  parameters  as  motor  diameter, 
launcher  length,  and  propellant  burning  rate  are  dealt  with. 
To  perform  an  optimized  design  study  OVERLAY  (l,i)  is  called 

while  OVERLAY  (1,2)  is  called  if  a certain  design  is  to  be 

i 

evaluated  with  a par«uneter  study.  OVERLAY  (1,2)  contains 

■i  ‘ 

an  advanced  graphics  package  for  plotting  all  results  if 
desired. 
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The  recoilless  rocket  system  model  is  contained 
in  OVERLAY  (2,0).  Such  parameters  as  recoil-mass  weight, 
launcher  length,  and  propellant  burning  rate  are  typical 
parameters  dealt  with  in  OVERLAY  (2,0).  To  conduct 
an  optimized  design  study  OVERLAY  (2,1)  must  be  called 
while  OVERLAY  (2,2)  must  be  called  to  evaluate  a particular 
design,  conduct  parameter  studies  or  plot  all  results. 
OVERLAY  (2,2)  also  contains  the  advanced  graphics  routines 
mentioned  above. 

The  flight  characterisitcs  of  both  the  conventional 
and  recoilless  syst'^ms  are  assumed  to  be  the  same  after  the 
missile  leaves  the  launch  tube.  This  allows  the  flight 
characteristics  of  the  missile  to  be  evaluated  in  OVERLAY 
(3,0).  The  aerodynamic  force  and  moment  coefficients, 
guidance  system  gain  setting  and  miss-distance  are  typical 
parameters  dealt  with  in  OVERLAY  (3,0).  As  in  the  other 
overlays,  OVERLAY  (3,1)  is  called  to  conduct  an  optimized 
design  study  while  OVERLAY  (3,2)  is  called  if  one  wishes  to 
conduct  a parameter  study  with  a graphics  capability. 

The  designer,  of  course,  does  not  need  to  understand 
the  overlay  system  nor  how  to  call  them  since  these  opera- 
tions are  conducted  automatically  as  a result  of  the  answers 
provided  to  certain  simple  questions.  These  characteristics 
are  illustrated  by  example  in  Section  7 of  this  report. 
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The  brief  descriptions  contained  in  this  section 
are  intended  to  mainly  provide  an  overview  of  the  CADLAW 
system.  Details  of  the  various  mathematical  models  used 

I 

to  describe  the  different  systems  are  contained  in  Sections 
3,  4,  and  5. 

2.2  Optimization  Strategy 

The  criterion  used  to  determine  what  constitutes  an 
"optimum  design"  must  be  decided  upon  by  the  designer.  The 
existing  CADLAW  system  has  certain  built-in  criteria  in  the 
form  of  specified  objective  functions.  However,  these  cri- 
teria may  be  changed  at  the  discretion  of  the  designer 
should  the  desire  arise. 

I 

There  are  two  types  of  objective  functions  employed 
in  CADLAW:  one  to  govern  the  optimization  of  the  propulsion 
phase  of  the  rocket  performance  and  one  to  regulate  the 
optimization  study  of  the  flight  phase.  The  optimum  design 
of  the  propulsion  phase  parameters  for  both  the  conventional 
and  recoilless  systems  are  based  upon  an  objective  function 
of  the  form 

''I'W  (2.1) 

r 

where  the  notation  Fj^(V^Wg)  indicates  the  objective  function 

i 

is  a function  of  the  muzzle  velocity  ' ^ and  system  weight 
Wg.  The  muzzle  velocity  and  system  weight  Wg.  The 
quantities  C^,  C2>  a,  and  3 are  constant  parauneters  to  be 
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specified  by  the  designer.  The  objective  function  in 
Eq.  (2.1)  was  selected  on  the  basis  of  convenience  and  what 
appears  to  be  a rational  design  philosophy.  The  importance 
of  muzzle  velocity  relative  to  system  weight  can  be  changed 
by  merely  changing  the  quantities  C^,  C2/  01,  and  B.  An 
optimum  design  is  defined  as  the  set  of  design  variables  which 
produce  a minimum  value  for  the  objective  function  F^.  This 
results  in  a system  with  low  weight  and  high  muzzle  velocity. 

The  design  variables  associated  with  the  flight  phase 
of  the  rocket  system  performance  are  selected  on  the  basis 
of  minimizing  the  miss  distance  d^  and  the  time  of  flight 
t^.  The  objective  function  selected  to  accomplish  this  is 
of  the  form 

F,(t.,d  ) = C.[t.]^  + C.[d  (2.2) 

where,  as  before,  C^,  C^,  A,  and  p are  parameters  with  which 
to  change  the  relative  importance  of  the  time  of  flight  and 
miss  distance. 

CADLAW  contains  default  values  for  the  parameters 
C2»  C^,  C^,  a,  3/  A,  and  n if  the  designer  does  not  wish  to 
specify  these  parameters  himself.  The  optimization  strategy 
is  to  search  combinations  of  design  variables  which  result 
in  a maximum  or  minimum  value  for  the  objective  function. 
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2.3  Design* Variables  and  Parameters 

One  concept  which  is  basic  to  the  optimization 
capability  of  CADLAW  concerns  design  variables  and  design 
parameters.  Obviously , there  are  literally  hundreds  of 
quantities  to  be  specified  in  calculating  the  propulsion 
and  flight  characteristics  of  a light  antitank  missile  sys- 
tem. State-of-the-art  optimization  methods  do  not  permit 
the  evaluation  of  a large  number  of  parameters  in  ah  optim- 
ized design  process.  Consequently,  the  design  quantities 
are  divided  into  two  groups,  design  variables  which  are 
selected  in  the  optimization  process,  and  design  parameters 
which  must  be  arbitrarily  selected  by  the  designer.  Al- 
though the  design  parameters  are  not  optimized  automatically, 
they  may  be  changed  in  such  a way  that  their  influence  can 
be  exeunined  and  appropriate  values  selected  for  the  design. 


3.  CONVENTIONAL  ROCKET  SYSTEM 

The  computer-aided  preliminary  design  system 
applicable  to  conventional  rocket  launcher  systems  assumes 
the  system  to  consist  of  a small  rocket  propelled  missile 
launched  from  a tube.  Figure  3.1  contains  a schematic 
diagram  of  the  system.  As  can  be  seen,  the  rocket  motor 
consists  of  a combustion  chamber  where  solid  propellant  is 
burned  and  a nozzle.  The  burned  gases  are  exhausted  from 
the  rear  of  the  laur.'ih  tube.  For  protection  of  the  person 
holding  the  launcher,  combustion  is  completed  and  the  cham- 
ber is  allowed  to  exhaust  to  a low  pressure  before  the  mis- 
sile leaves  the  launch  tube. 

Typical  parameters  of  interest  in  the  preliminary 
design  of  this  conventional-in-tube-burning  rocket  system 
are  the  system  diameter,  launcher  length,  propellant  burning 

I 

rate,  propellant  mass,  chamber  pressure  and  nozzle  expansion 
ratio.  CADLAW  allows  these  design  variables  to  be  selected 
such  that  the  system  performance  is  optimized. 

As  discussed  in  Section  2.3,  the  performance  of  the 
conventional  rocket  system  is  measured  by  the  value  of  an 
arbitrarily  defined  objective  function  which  depends  on  the 
muzzle  velocity  and  system  weight.  An  optimum  design  is 
defined  as  a set  of  design  variable  values  which  minimizes 
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the  objective” function  and  thereby  maximizes  the  muzzle 
velocity  while  keeping  the  system  weight  low.  To  calculate 
the  muzzle  velocity  and  system  weight  when  all  design  para- 
meters are  specified,  resort  must  be  made  to  a system  of 
mathematical  relations  obtained  by  application  of  the 
principles  of  mechanics. 

3.1  Ballistic  Performance  Equations 

The  mathematical  equations  which  describe  the  ballistic 
performance  of  small  solid  rockets  are  described  in  detail  in 
the  literature  [14] . For  the  sake  of  completeness,  the 

equations  are  repeated  here  in  the  form  and  order  in  which 
they  are  evaluated  in  CADLAW. 

The  following  parameters  are  assumed  to  be  specified 
numerically  at  the  time  the  performance  equations  are  evalu- 
ated to  obtain  the  muzzle  velocity  and  system  weight: 

1.  = chamber  pressure  at  ambient  temperature 

(psi) 

2.  A^/A^  = nozzle  expansion  ratio  (dimensionless) 

3.  = outer  diameter  of  motor  case  (in.) 
m 

4.  kAj^  = propellant  volumetric  burning  rate 

3 

coefficient  (in  /sec) 

5.  Wp  = weight  of  solid  propellant  (lbs) 

6.  L = length  of  motor  combustion  chamber  (in.) 

m 

7.  a = specific  heac  of  gas  (BTU/lb“F) 

3 

8.  p = propellant  density  (Ibs/in.  ) 
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9. 


10. 

11. 

12. 

13. 

14. 

15. 

16. 


FS 

T_ 


IT, 


m 


17. 

18. 

19. 

20. 


a 


w. 


w 


W. 


system  factor  of  safety 
cold  temperature  extreme  (®P) 
ambient  temperature  (®F) 
hot  temperature  extreme  (®F) 
pressure  sensitivity  parameter  (%/®F) 
final  chamber  blowdown  pressure  (psi) 

thickness  of  launcher  (in.) 

density  of  motor- launcher  material 
3. 


(Ibs/in.'') 

= allowable  stress  for  motor  material  (psi) 
= nozzle  expansion  angle  (degrees) 

= warhead  weight  (lbs)* 

- weight  of  fins,  pole  piece,  fixtures, 

Wik 

etc.  (lbs) 

With  numerical  values  specified  for  the  above  parameters,  the 
following  geometrical  quantities  can  be  evaluated: 


H = (5;)  (Ae> 


and 


c 3 


(3.1) 

(3.2) 


(3.3) 
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where  A is  the  area  of  the  nozzle  exit  plane,  A,  is  the 
e t 

area  of  the  nozzle  throat  and  V is  the  volume  of  the  motor 

c 

combustion  chamber. 

The  Mach  number  at  the  nozzle  exit  plane  can  be 
evaluated  by  solving  the  nonlinear  equation 


2 

(Y-1) 


Y+1 

2 


<«e^> 


2 (Y-1) 
(Y+1) 


1] 


(3.4) 


for  Mg.  This  is  accomplished  in  CADLAW  by  using  a rapid 
step-by-st  evaluation  routine.  Based  upon  a knowledge  of 
the  motor  chamber  pressure  P at  ambient  temperature  T,  and 
the  temperature  extremes  and  T^,  the  chamber  pressures  at 
the  temperature  extremes,  and  ^ch  , can  be  evaluated  as 


P 


cc 


P 

ca 


200  + TT,  (T  -T  ) 
k c a 

200  -•  IT,  (T  -T  ) 
k c a 


(3.5) 


and 


P 

ca 


200  - 


(3.6) 


The  conventional  system  performance  varies  with  temperature. 
For  a given  set  of  parameters  the  system  muzzle  velocity  and 
chamber  pressure  will  be  lower  when  operated  at  the  cold  tem- 
perature extreme  than  when  operated  at  the  hot  temperature 

! 

extreme.  The  system  weight  depends  strongly  on  the  motor 
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chamber  pressure.  For  conservative  purposes  the  system 
muzzle  velocity  and  launcher  weight  are  based  upon  opera- 
tion at  the  cold  temperature  extreme  while  the  motor 
weight  is  based  upon  the  chanber  pressure  when  the  system 
is  operated  at  the  hot  temperature  extreme 

The  pressure  at  the  nozzle  exit  plane  for  tem- 

perature is 


P 

ec 


cc 


1 + 


Y-1 

2 


Y 

1-Y 


(3.7) 


The  thrust  F at  T is  given  by  the  equation 
c c * 


F = 
c 


Iti  2 

Y-1  2Y 
Y-l' 


[ 1 


1/2 

X 


[1/2  (1  + cosa) ] 


(3.8) 


and  the  propellant  mass  discharge  rate  at  temperature 

T by  the  relation 
c 


W = p kA,  (P  )^  = 
c . p b cc' 


P Kw(P  ) 
^p  b cc 


(3.9) 


where  the  parameter  one  of  the  optimum  design  variables 

in  CADLAW,,  Knowing  the  thrust  and  mass  discharge  rate  per- 
mits the  calculation  of  an  effective  exhaust  velocity  V by 
use  of  the  equation 
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V. 


ec 


(3.10) 


w^/g 


The  burning  time  of  the  propellant  tj^  can  be  determined  by 
the  equation 


t = ^ 
be 


(3.11) 


where  Wp  is  the  weight  of  the  propellant.  Using  the  laquation 


W 


®d  A^P 


(3.12) 


t cc 


to  calculate  the  mass  discharge  coefficient  permits  the 
elapsed  time  from  the  time  when  burning  is  completed  until 
the  motor  exhausts  to  a pressure  P^,  i.e.,  the  "blowdown" 
time  tj^^,  to  be  calculated  by  using  the  equation 


V C , 
c d 


'bd  6gY(Y“l)A^ 


Yfl  ^ 
(^)  Y-1 


(^)  2y  - 1 

L 


Y-1 


(3.13) 


The  chamber  pressure,  thrust,  propellant  burning  time, 
and  motor  blowdown  time  have  been  calculated  using  flqs. 

(3.1 — 3.13).  Before  the  missile  exit  velocity  can  be  deter- 
mined, the  total  missile  weight,  both  with  propellant  and 
without  propellant  W , must  be  known.  This  includes  the 
structural  weights  of  the  motor  and  nozzle.  A description 
of  the  structural  and  system  weight  calculations  are  contained 
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in  Section  3.2.  It  is  sufficient  to  recognize  that,  from 

a performance  viewpoint,  the  structural  weight  depends 

directly  on  the  highest  operational  pressures  '■oh-  The  other 

weight  components  are  specified  by  the  designer. 

Assuming  the  weights  W , and  W are  known,  the  mis- 

wi  wu 

sile  velocity  at  the  time  of  burnout  v^^  is 


V = V In  (W„,/W  ) 

m e wi  wu 


(3.14) 


whjsre  ln(x)  denotes  the  natural  logarithm  of  x.  The  specific 
impulse  can  be  expressed  as 


I 


sc 


(3.15) 


and  the  distance  traveled  during  burning  Sj^  as 

The  distance  traveled  during  blowdown  can  be  expressed  as 


^bd  “ '^m  * ^bd 


(3.17) 


since  the  missile  velocity  is  assumed  to  reamin  constant 
from  the  time  of  burnout  until  the  missile  exits  the  launch 
tube.  Therefore,  v^^^  is  the  muzzle  velocity. 
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3.2 


Structural  Design  and  Weight  Equations 


As  indicated  in  Section  3.1,  the  structural  design 


and  resulting  weight  depends  directly  on  the  maximum  opera- 


ting pressures.  The  effects  of  high  acceleration  are  small. 


relative  to  the  pressure  induced  stresses  and  strains,  due 


to  the  small  size  of  the  weapon  systera.  The  design  pressure 


is  obtained  by  multiplying  the  maximum  chamber  pressure 


by  thf3  factor  of  safety  FS,  i.e.. 


"■d  ' ^oh 


(3.18) 


Thin-walled  pressure  vessel  theory  is  utilized  in  the 


structural  design  oa  the  rocket  motor  case  and  nozzle.  The 


rocket  motor  structure  is  idealized  as  illustrated  schemati- 


cally in  Figure  3.2.  The  nozzle  consists  of  two  truncated 


cones  of  thickness  t^,  the  chamber  is  a thin-walled  circular 


cylinder  and  the  dome  consists  of  a 2:1  ellipse  of  thickness 


Based  upon  a maximum  normal  stress  failure  criterion. 


the  weight  of  the  nozzle  can  be  shown  to  be 


^'n  = ^ 


3 2 

[4(^)  + - 7 .5833] 


(3.19) 


where  c is  the  maximum  allowable  stress  in  the  nozzle.  The 
a 


cylindrical  chamber  weight  can  likewise  be  expressed  as 
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.WWI 


Trp  L D ‘ P, 

rciniin  r ^ i 

2 


(3.20) 


and  the  end-closure  dome  weight  expressed  as 


'^d  " 128  -P-  ' ^ ^ °a~^d 

4 d a a 


(3.21) 


the  unloaded  motor  structural  weight  can  be  expressed  as 


W = W + W +W,  + W (3.22) 

m n c d X 

where  is  the  weight  of  the  non-structural  components  such 
as  the  polar  adaptor  and  guidance  fins. 

The  unloaded  weapon  weight  is  calculated  by  adding  the 
weight  of  the  warhead  to  the  unloaded  motor  weight,  i.e.. 


W 


wu 


= + 
m 


w 


w 


(3.23) 


The  initial,  loaded  weapon  weight  is  obtained  by  adding 
the  weight  of  the  propellant  Wp  to  the  unloaded  weapon  weight 


W. 


wl 


= W + 
wu 


W 


(3.24) 


In  addition  to  the  motor  structural  weight  and  the 
warhead  weight,  the  system  weight  must  include  the  weight  of 
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the  launcher.*  The  launcher  length  depends  upon  the  burning 
time  of  the  motor  and  the  distance  traveled  during  burning 
and  subsequent  blowdown.  These  distances  have  been  calculated 
in  Eqs.  3.14  and  3.15.  With  numerical  values  for  Sj^  and 
the  launcher  length  can  be  expressed  as 

Lj  = (3.25) 

and  the  launcher  weight  as 

= "ViLi.  1 (3-26) 

where  t^^  is  the  thickness  of  the  launch  tube  as  specified 
by  the  designer.  Finally,  the  total  system  weight  W is  the 
sum  of  the  loaded  weapon  weight  and  launcher  weight 

= W , + W,  (3.27) 

s . wl  1 

As  discussed  earlier,  the  purpose  of  the  mathematical 
modeling  of  the  conventional  rocket  performance  is  to  permit 
the  muzzle  velocity  v^^  and  system  to  be  calculated  as  a 
function  of  the  design  variables,  kj^,  d^^^,  L^,  and  the 
design  parameters.  The  muzzle  velocity  and  system  weight 
have  been  calculated  in  Eqs.  2.14  and  2.27,  respectively. 

Their  dependence  on  the  design  variables  can  be  seen  to  be 
quite  complicated.  This  complication  precludes  the  use  of 
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■ 


-•1  ■T7»- 


any  analytical,  optimization  technique  and  instead  requires 
the  use  of  a pattern  search  routine  such  as  that  described 
in  Section  6. 


4.  RECOILLESS  LAUNCHER  SYSTEM 


A schematic  diagram  of  the  recoilless  launcher 
system,  upon  which  the  mathematical  model  is  based,  is  il- 
lustrated in  Figure  4.1.  The  parameters  utilized  in  the 
analysis  of  the  recoilless  launcher  system,  some  of  which 
are  contained  in  Figure  4.1,  are  defined  as  follows: 

= weight  of  warhead,  (lbs) 

W2  = weight  of  recoil  mass,  (lbs) 

L = total iilength  of  launch  tube,  (in.) 

* portion  of  launch  tube  traversed  by  warhead, 

(in. ) 

L2  = portion  of  launch  tube  traversed  by  recoil  mass, 
(in. ) 

W = weight  of  solid  propellant  charge,  (lbs) 

hr 

d = inside  diameter  of  launch  tube,  (in.) 
t^  = wall-thickness  of  launch  tube,  (in.) 
tj.  = thickness  of  reinforcing  ring,  (in.) 

L^  = length  of  reinforcing  ring,  (in.) 

K = .propellant  burning  rate  coefficient, 

3 

(in.  /Ib-sec/sec) 

n - burning  rate  exponent  (dimensionless) 

= initial  position  of  warhead,  (in.) 
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^2  ^ initial  position  of  recoil  mass,  (in.) 

The  objective  of  the  analysis  is  to  develop  a means 
for  calculating  the  exit  velocity  of  the  warhead  and  the 
total  system  weight  when  all  parameters  of  the  system  are 
defined.  The  performance  of  the  system  can  then  be  evaluated 
in  terms  of  the  warhead  exit  velocity  and  system  weight 
values.  An  optimum  design  is  defined  as  one  which  simultan- 
eously maximizes  the  muzzle  velocity  and  minimizes  the  system 
weight,  according  to  some  predetermined  relation. 

To  calculate  the  muzzle  velocity  and  system  weight, 
resort  must  be  made  to  various  principles  of  mechanics. 
Specifically,  the  principles  of  balance  of  momentum,  mass, 
and  energy  must  be  employed  along  with  concepts  of  structural 
mechanics  and  fluid  dynamics  to  arrive  at  mathematical  equa- 
tions which  describe  the  system  behavior. 

4.1  Equations  of  Motion 

Application  of  Newton's  Second  Law  yields  the  equa- 
tions of  motion  for  the  warhead  and  recoil  mass: 


and 


(4.1a) 


(4.1b) 
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where 

= weight  of  warhead,  (lbs) 

W2  = weight  of  recoil  mass,  (lbs) 

= position  of  warhead  at  time  t,  (in.) 

X2  = position  of  recoil  mass  at  time  t,  (in.) 

Pg  = pressure  of  gas  in  chamber,  (psi) 
d = inside  diameter  of  launch  tube,  (in.),  and 
g = gravitational  constant,  (in. /sec  ) . 

A dot  over  a variable  is  used  to  denote  the  derivative  of 

th.at  variable  with  respect  to  time.  The  initial  conditions  j 

are  i 

I 

• ( 

I 

1)  @ t = 0,  = A,  X2  = 0 and  (4.2a)  . . 

2)  0 t = 0,  Xj^  = X2  = 0 (4.2b) 


where  A is  the  initial  position  of  the  warhead.  From  Eqs. 
(4.1a)  and  (4.1b),  it  is  observed  that  the  acceleration  of 
the  warhead  and  recoil  mass  are  related  as 


Integration  of  Eq.  (4.3)  and  employment  of  the  initial  con- 
ditions listed  in  Eqs.  (4.2a)  and  (4.2b)  allows  the  establish- 
ment of  a linear  relationship  between  the  respective  positions 
of  the  two  masses  and  between  the  forward  and  rearward  lengths 
of  the  launch  tube.  These  relationships  are 
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*2  *“  Wj 

and 

h " ^ 'V'^'  '^•=> 

where  A Is  the  initial  position  of  the  warhead  with  respect 
to  the  origin  of  the  coordinate  system  illustrated  in  Figure 

4.1. 

The  relationship  expressed  in  Eg.  (4.4)  eliminates  the 
necessity  of  solvinc  Eg.  (4.1b).  By  solving  only  the  eguatlon 
of  motion  for  the  warhead,  Eg.  (4.1a),  the  motion  of  the  re- 
coil mass  can  be  determined  through  the  use  of  Eg.  (4.4). 

As  will  be  shown  in  subseguent  developments,  the  eguation  of 
motion  for  the  warhead  is  highly  nonlinear  and  must  be  solved 
using  numerial  integration  procedures.  To  this  end.  Eg.  (4.1a) 
can  be  replaced  by  the  following  two  first  order  differential 
eguations : 


(4.6) 


and 


dv, 

dt  “ 4Wj^ 


(447) 
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where  is  the  speed  of  the  warhead  in  units  of  in. /sec. 

4.2  Continuity  Equation 

The  principle  of  conservation  of  mass  can  be  expressed 
here  by  equating  the  mass  rate  of  propellant  combustion  to 
the  mass  rate  of  gas  generation  within  the  pressure  chamber, 

i.e. , 


= ^(«gl  = 


(4.8) 


where 

Wg  = gas  weight,  (lbs) 

r = propellant  burning  rate,  (in. /sec.) 

2 

= propellant  burning  surface  area,  (in.  ) 

3 

p = propellant  density,  (Ibs/in.  ) 

Jr 

p = gas  density,  (Ibs/in.', 

y 

and 

3 

V = volume  of  gas  in  chamber,  (in.  ) 

9 

The  volume  occupied  by  the  gas  at  any  time  can  be  expressed 
in  terms  of  the  distance  between  the  warhead  and  the  recoil 
mass  and  the  launch  tube  diameter  or,  using  Eg.  (4.4), 


+ X2) 


(4.9) 
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or,  using  Eq.  (4.4)  , 


„ i.a^  ,,V«2,  «i„ 

''g  ' — ' ‘-Wj-’  '‘l  - 
Thus,  the  weight  of  the  gas  at  any  time  is 


(4.10) 


_-2  W,+W_  W. 

"g  “ ‘’g''g  = V OgO^’^V  wT^> 


(4.11) 


These  equations  must  be  further  developed  through  the  use  of 
the  propellant  burning  rate  and  gas  constitutive  equations. 

4 . 3 Propellant  Burning  Rate 

and  Gas  Constitutive  Equations 

The  propellant  burning  rate  is  assumed  to  depend  on  the 
chamber  pressure  through  the  relation 

r = kPg”  (4.12) 

where  k and  n are  constants  and  the  gas  is  assumed  to  obey 
the  perfect  gas  law,  i.e., 

Pg  . P^RT^  (4.13) 

In  Eq.  (4.11),  R is  a gas  constant  (37.69908  in./®R)  and  T^ 
is  the  local  gas  temperature,  in  degrees  Rankine.  Furthermore, 
the  tempersture-pressure  relationship  assumed  for  this  gas  is 


2T  , 

='^-tan"  (P  /P„) 
g TT  ' g'  o 


(4.14) 


Where  T^  and  P^  are  reference  temperature  and  pressure,  re** 
spectively.  T^  is  assumed  to  have  the  value  of  5500®R  and 
the  reference  pressure  P is  100  psi. 


33 


% 


4. 4 Combined  Field  Equations 

Through  the  use  of  Eqs.  (4.13)  and  (4.14),  the  gas 
weight  can  be  expressed  in  terms  of  the  pressure  and 

I 

position  as 


, g 

^ 8RT^  W2  tan’^(Pg/P^) 

Wf 

“ W^+W2^^ 


(4.15) 


For  simplicity  of  notation,  if  we  define  two  constants  a and 
0 as 

Wji 


and 

8RT  W2P  A k (4.16) 

P ■ 2 i — 

■n  d (Wj^+W2) 


then  the  conservation  of  mass  equation,  Eq.  (3.6)  can  be 
expressed 


(Xj^-a) 


tan 


-1 


3 1 


(4.17) 


Expanding  the  derivative  in  this  equation  gives 
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3P 


« v,P_ 

n 1 g 


9 tat.-l  (P^/P^) 


+ (x^-a)Pg  X 


1 P /P 

1 _ q O 

tan"^(P  /P^)  11+  (P„/P«) [tan’^ (P«/P«)  1 ^ 


(4.18) 


g'  O' 


g'  o 


g'  o' 


where,  as  in  Eqs.  (4.6)  and  (4.7),  is  the  speed  of  the 
warhead.  Eq.  (4.18)  represents  a nonlinear  first  order  dif- 
ferential equation  in  terms  of  the  pressure  P^. 

Solving  Eq.  (4.18)  for  P^  and  restating  Eqs.  (4.6)  and 
(4.7)  allows  the  final  governing  equations  to  be  summarized 
as  three  simultaneous,  first  order,  nonlinear  differential 
equations  of  the  form 


Xi  = v^ 


(4.6) 


V = p 

'^1  4Wj^  ^g 


(4.7) 


and 


3P, 


n 


(Xj^-a) 


tan-l  (Pq/P^) 

(pITpIT" 


(4.19) 


g o 

tan"^(P  /P^)  (1+ (P  /P^) [tan"^ (P„/P^) 1 ^ 


g'  o' 


g'  O' 


g'  o 


I 


In  addition,  the  rate  of  propellant  combustion  is  determined 
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by  the  differential  equation 

iL  = - )cA.  P ” (4.20) 

P b p g 

This  equation  must  be  solved  to  determine  the  time  at  which 
propellant  burning  ceases.  The  initial  conditions  which 
apply  for  the  above  equations  are: 


and 


1)  = A 

2)  v^  = 0 


3) 


g o 


4) 


p o 


(4.21) 


4 . 5 Conditions  After  Propellant  Burnout 

The  governing  Eqs.  (4.6),  (4.7)  and  (4.19)  apply  only 
until  the  propellant  is  totally  consumed.  After  this  time 
of  burnout,  the  gas  in  the  pressure  chamber  is  assumed  to 
expand  isentropically  until  the  warhead  and  recoil  masses 
leave  the  launch  tube.  The  equations  of  motion  of  the  warhead 
and  recoil  masses  are  the  same  as  before  burnout,  Eqs.  (4.1a) 
and  (4.1b),  and  the  positions  of  the  two  masses  can  be 
related  as  in  Eq.  4.4.  For  isentropic  expansion  the 
pressure  depends  on  the  temperature  as 


Y 


k,T 

3 q 


Y-1 
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(4.22) 


where  Is  a constant  to  be  determined  from  the  cohditlons 
at  burnout  and  y is  the  ratio  of  specific  heats.  After  burn- 
out the  weight  of  the  gas  is  constant 


TTd' 


W. 


(Wj^+W2) 


Y-1 


(k^) 


1 

Y 


(Xj^-a) 


(4.23) 


where  all  of  these  parameters  have  been  defined  earlier.  By 

1/y 

inspection  of  Eq.  (4.23)  the  product  ''  (Xj^-a)  is  assiuned 
to  be  a constant 


(Xj^-a) 


lA  * 

(Xj^  -a) 


(4.24) 


n ★ 

Where  and  are  the  gas  pressure  and  position  of  the 
warhead  at  the  time  of  burnout.  Thus,  P^  can  be  expressed 
as 

<4.25) 

g g 1 I 


and  the  equation  of  motion  becomes 


(^)P_  = (^)  P*  (x^*-a)^  (x^-a)-^  = !^xj_ 


(4.26) 


This  equation  can  be  integrated  in  closed  form  if  it  is  first 
expressed  in  terms  of  the  velocity  as 
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dx 


2 ^ 
Tid‘^gP„ 

g 

4W, 


* Y 

(x,  -ar 


(Xj^  - a) 


“Y 


(4.27) 


Integrating  each  side  of  this  equation  and  denoting  the 
warhead  velocity  as  when  x=Lj^/  allows  the  exit  velocity  to 
be  expressed  as 


2 * * 
gTrd  Pg  (x^  -a) 

2Wj^  (Y-1) 


Y 


1 - 


■(x^  -a)! 

(L^-a) 


,Y-1- 


(4.28) 


it 

where  it  is  understood  that  x^^>  a and  a. 

' » 

4. 6 Solution  Method 

As  mentioned  earlier,  the  object  of  this  development 
is  to  permit  calculation  of  the  warhead  and  recoil  masses 
exit  velocities  and  the  system  weight.  The  exit  velocities 
of  the  two  masses  can  be  calculated  from  Eqs.  (4.28)  and 
(4.4)  and  the  system  weight  can  be  calculated  from  a know- 
ledge of  the  peak  pressure.  To  evaluate  Eq.  (4.28)  and  the 
peak  pressure,  the  governing  field  equations,  Eqs.  (4.6), 
(4.7),  and  (4.19)  must  be  integrated  with  respect  to  time. 

A fourth  order  Runge-Kutta  integration  procedure  has 
been  utilized  to  solve  this  sytem  of  equations.  To  allow 
for  a description  of  this  solution  method>  let  the  govern- 
ing equation  be  expressed  in  vector  form  as 


i 


d 

dt 


X = F(t,x) 


(4.29) 
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T 

where  x = [x, , v , P ] and  F is  the  right  hand  of  the 

X X y •• 

governing  equations  and  may  depend  on  Xj^,  or  P^.  By 
discretizing  the  field  variables  to  correspond  to  discrete 
time  values  t^^,  separated  by  a time  increment  At,  the  Runge- 
Kutta  method  provides  an  accurate  method  for  calculating  the 
vector  X which  contains  the  field  variables  at  bime 
t^^j^  if  the  values  of  x^  at  time  t^  are  known.  The  algorithm 
can  be  expressed  as 

x^'*'^  = x^  + — [Ki  + 2K2  + 2K3  + K^]  , (4.30) 

where 

K,  = F(t.,x^)  (4.31) 


and 


52  = !<*i  + ¥'  t 5i> 
5a  = !<*i  * f ' 5'  ¥ 52) 

K4  = + ^^^3). 


(4.32) 

(4.33) 

(4.34) 


This  procedure  may  be  employed  to  determine  the  values 
* * ■ * 

of  Xj^,  Vj^,  and  P^  at  the  time  of  burnout  as  well  as  the  maxi- 
mum pressure  in  the  cheimber.  With  these 'values  known  the 
velocities 'of  the  two  masses  can  be  determined  from  Eqs.  (4.28) 
and  (4.4)  and  the  system  weight  determined  from  the  peak  pres- 
sure calculation. 
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4.7  System  -Weight 

The  system  weight  W may  be  expressed  in  terms  of 

s 

system  component  weights  as 

”s  = "l  ”2  + It  + ”p  + ”x 

where  and  V?2  are  the  weights  of  the  warhead  and  recoil 
masses,  is  the  weight  of  the  launch  tube,  is  the  pro- 
pellant weight,  and  is  the  weight  of  such  extras  as 
handles,  straps,  sights,  etc.  The  weight  of  the  launch  tube 
can  be  calculated  assuming  a maximum  normal  stress  failure 
criterion  and  using  thin-walled  tube  theory  to  calculate  the 
maximum  stress  in  the  tube  as  a function  of  the  maximum  in- 
ternal pressure  P . P must  be  determined  by  integrating 

max  max 

the  nonlinear  differential  equation,  Eq.  (3.19),  as  dis- 
cussed earlier.  The  maximum  pressure  is  selected  from  the 
entire  pressure  history  during  burning.  A design  pressure 
for  the  launch  tube  is  computed  as 

'■a  = ‘■max 

Design  of  the  launch  tube  requires  a special  feature 
near  the  initial  position  of  the  warhead  and  recoil  mass. 

Very  early  after  the  propellant  is  ignited  and  before  the  war- 
head and  recoil  mass  have  had  a chance  to  move  significantly, 
the  burning  gas  pressure  within  the  launch  tube  increases  to 
an  extremely  high  pressure.  After  the  two  masses  begin  to 
move,  the  pressure  within  the  launch  tube  drops  considerably. 
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If  the  weight *o£  the  launch  tube  were  designed  based  upon 
these  early  high  pressures,  the  system  weight  would  be  un- 
acceptably high.  To  avoid  this  problem  and  thus  reduce  the 
system  weight,  the  launch  txibe  is  assumed  to  have  a high- 
strength  reinforcing  ring  which  surrounds  the  launch  tube 
in  the  vicinity  of  the  initial  position  of  the  masses.  This 
feature  is  illustrated  in  Figure  4.1.  The  reinforcing  ring 
weight  can  be  expressed  as 


W 


ird^P,  (L  )p 
L d r r 


(4.37) 


where  L^.  is  the  length  of  the  ring,  normally  assumed  to  be 

* 

approximately  3 inches,  and  is  the  ultimate  strength  of 
the  reinforcing  ring. 

Finally,  the  weight  of  the  launch  tube  can  be  expres- 
sed as 


Wj,  = «r  * 


ndL  P^(Lj+L2)Pj_ 


a 


(4.38) 


where  a is  the  maximum  stress  allowable  in  the  non-reinforced 

cl 

portion  of  the  launch  tube.*  Eq^  (4.38)  now  permits  the  entire 
system  weight  to  be  calculated  as  in  Eq.  (4.35). 

The  muzzle  velocity  v^  and  system  weight  can  be 
calculated  as  indicated  in  Eqs.  (4.28)  and  (4.38),  respec- 
tively. Obviously  it  is  not  possible  to  express  the  dependence 
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of  these  quantities  on  the  design  variables  in  any  sort  of 
functional  form.  This  precludes  the  use  of  an  analytical 
optimization  routine  and  instead  requires  that  a pattern 
search  method  such  as  that  described  in  Section  6 be 
utilized  to  select  an  optimum  design. 
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5.  TRAJECTORY  ANALYSIS 

To  determine  the  effectiveness  of  light  antitank 
weapons,  the  capability  for  hitting  the  target  must  be 
evaluated  by  computing  the  trajectories  of  the  missile  and 
target.  This  calculation  involves  modeling  the  aerodynamic 
and  ballistic  forces  on  the  missile  as  well  as  any  guidance 
and  control  systems  which  exist  and  integrating  the^e  sys- 
tems into  a mathematical  model.  Due  to  the  exploratory 
nature  of  this  projtict,  a trajectory  analysis  was  developed 
which  embodies  all  the  pertinent  featureis  of  the  system 
but  allows  for  guidance  and  control  characteristics  only 
within  a horizontal  plane.  Motion  normal  to  this  plane  was 
assumed  to  follow  a ballistic  trajectory.  This  assumed  motion 

I 

corresponds  closely  to  the  motion  of  light  antitank  weapons 
due  to  the  high  muzzle  velocities  and  to  the  character  of 
the  evasive  actions  of  the  targec.  The  trajectory  analysis 
includes  target  motion  and  random  crosswind  effects. 


5.1  Equations  of  Motion 

Th6  coordinate  system  and  system  parameters  used  to 
define  motion  of  the  missile  are  illustrated  in  Figure  5.1. 
Control  of  -the  missile  is  maintained  with  a rudder.  The 
equations  of  motion  of  the  missile  can  be  expressed  as 
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(5.1) 


and 


ZF 


= -m  V 


ZM  = 


(5.2) 


where 

F = force  vector,  (lbs) 

M = moment,  (in. -lbs) 

V = velocity  vector,  (in. /sec.) 

*1'  = angle  of  rotation,  (radians) 

2 

m = mass  of  missile,  (Ibs-sec  /in.) 

2 

I = mass  moment  of  inertia,  (in. -Ibs-sec  ) . 
zz 


The  velocity  vector  can  be  expressed  in  terms  of  the  unit 

vectors  e and  e through  the  following  equations 
X y 


V = (v)^  ^ Y' 


(5.3) 


W 


(5.4) 


and 


<Pvy) 


+ tv  )e 
X y 


(5.5) 


The  equations  of  motic  :an  then  be  expressed  in  component 


form  as 


EF  = m(v  - 4'v  ) 
X X y 


(5.6) 


EF  = m(v  + 4'v  ) 
y y X 


(5.7) 


and 

EM  = I 'i'  (5.8) 

z zz 

To  make  solution  of  these  equations  possible,  linearize  the 
velocity  components  about  an  equilibrium  position  as 


v„  = U + Au 
X o 

(5.9) 

Vy  = Av 

(5.10) 

and 

T =4'  + A4» 

o 

(5.11) 

where 

Au,  Av  and  A4*  are  pertubations 

about  the  equilibrium 

state. 

The  equations  of  motion,  Eqs. 

(5.6)  — (5.8),  reduce  to 

EF  = mAu, 
X ' 

(5.12) 

EF  = m (Av  + U A4'  ) , 

(5.13) 

y o 

and 


SMz  = . (5.14) 
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5.2  i Aerodynamic  Forces 

Aerodynamic  forces  are  exerted  on  the  missile  due  to 
drag  and  pitching  moment  effects.  As  illustrated  in  Figure 
5.1,  the  relative  wind  acts  at  an  angle  3 with  respect  to 
the  liiissile  axis.  The  relative  wind  can  be  expressed  in 
vectpr  form  as 


V y = V - w 
~m/w 


(5.15) 


or,  in  terms  of  the  previously  defined  quantities 


‘I 

V . = (U  +Au)  e + Avd  - w d - w § 

~m/w  o x y XX  y y 


(5.16) 


Vw  ' '“o*  ■ "x™®  ''  • *y  ''>®x 


+ (Av  + w sin  H*  - w„  cos'l')e 
X y y 


(5.17) 


Then,  3 can  be  expressed  as 


3 = tan 


-1 


Av  + w sin  'F  - W cos  *F 
X X 


U + Au  - W cos  V *-  W„  sinf 
o X y 


(5.18) 


If  3 ‘and  'F  are  assumed  small 


- o 


(5.19) 


The  aerodynamic  forces  F , F and  M are  functions  of  the 

X Y Z 
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A 

dynamic  pressure  q,  a characteristic  area  S,  the  missile 
length  b,  as  well  as  the  variables  Au,  3»  6^.  and  V.  For 
convenience,  introduce  the  dimensionless  velocity  components 
Au  and  Av  defined  as 


Au  = U Au  and  Av  = U Av 
o o 


and  express  the  aerodynamic  forces  as 


b 


fy' 9S(Cyge  >' 


and 


= qSb(C^33  t A^t 


(5.20) 


(5.21) 

(5.22) 

(5.23) 


where 


<J  = 5 %/w 

q ^ Ip  [(Uq  “ W^)2  + Wy^],  (5.25) 


and  where  the  coefficients  are  force  and  moment  coefficients 
detemined  by  experiment  or  estimated  from  tests  on  similar 
vehicles. 

The  equations  of  motion  can  then  be  expressed  as  I 
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■rw)n?^ 


* ' 

qr  = <=x  + ‘^xu'^"  + =x$  » * W'  ' ‘®-“’ 

^ (Av  + At)  =0.6  + 2^  C A Af  + C a «^.  (5.27) 

qS  ^ o ^ 

and 


qlE  « = =ze  ^ ^ 


b 

2U. 


'z5^ 


(5.28) 


The  velocity  components  of  the  missile  relative  to  the  (X,Y) 

<Cl 

reference  coordinate  system  can  be  expressed  as 


% 


U^[(l+Au)  cos  f - Av  sin  4'], 


(5.29) 


and 

•« 

Y = U [ (1+Au)  sin  f + Av  cos  V]  (5.30) 

m o 

These  governing  equations  must  be  integrated  in  time 
subject  to  the  initial  conditions 


X (o)  = 0,  Y„(0)  = 0,  Au(0)  = aG  , 
mm  o ^ 

Av(0j  - Av^r  A4'(0)  = V(0)  = W^. 


(5.31) 
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5.3  Solution  Technique 

The  differential  equations  of  motion,  Eqs,  (4.26) — 
(4.28),  represent  e set  of  coupled  second  order  linear  dif- 
ferential equations  which  can  be  reduced  to  a set  of  first 
order  differential  equations  of  the  form 

I^(aG)  = (5.32) 


(5.33) 


(5.34) 


dt 


(AV) 


(5.35) 


where  F = Av  + AT  and  the  other  quantities  in  Eqs.  (5.32) — 

(5.35)  can  be  inferred  by  comparison  with  Eqs.  (5.26) — (5.28). 
These  equations,  Eqs.  (5.32)  — (5.35),  can  be  cast  into  vector 
form  as 


where 

= [Au,  r.  Si,  AT] 

and 

T 

F = [g^^,  g2»  93/  • 


(5.36) 


50 


-'K\ 


This  vector  differential  equation  is  now  in  the  seune  form 
as  Eq.  (4.29)  and  the  step-by-step  Runge-Kutta  Integration 
method  described  in  Section  4 can  be  used  to  solve  Eq.  (5.36). 

With  the  solution  of  Eqs.  (5.32)  — (5.35)  known  at  each 
Instant  of  time,  Eqs.  (5.29)  and  (5.30)  can  be  solved  to  ob- 
tain the  X and  Y coordinates  of  the  missile  in  a step-by- 
step  process. 

5.4  Control  Law 

The  contol  law  used  in  CADLAW  to  Insure  intercept  of 
the  missile  with  the  moving  target  is  a simple  rate  dependent 
relation  of  the  form 

+ G2C  ■ (5.37) 


where 


? = y-  - Y, 


m 


(5.38) 


and 


5 


(5.39) 


and  where  y.  and  Y are  coordinates  of  the  target  and  missile, 
t m 

respectively..  The  parameters  and  G2  are  gain  settings 
and  can  be  determined  in  an  optimization  process. 


5.5  Target  Motion 

To  evaluate  the  performance  of  the  aerodynaunic  and 
guidance  and  control  systems,  the  antitank  missile  Is  assumed 
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to  be  attempting  to  intercept  a target  which  moves  in  the 
horizontal  plane.  To  provide  a variety  of  evasion  tactics, 
the  target  is  assumed  to  follow  a curve  defined  as 

2 2 

ax^  + + cy^  + dx^  + ey^  + f = 0 (5.40) 


where  (x^,y^)  denote  the  target  location  and  where  a,b,c  . . ., 
f are  arbitrary  constants  used  to  describe  the  target  notion. 
The  target  speed  v^  is  assumed  constant  and  may  also  be 
changed  by  the  designer  to  study  the  effects  of  different 
evasive  tactics.  Wi»th  the  target  speed  constant  the  target 
location  can  be  determined  at  any  time  hy  integrating  the 
differential  equations 

Xt  = V.  cos  6 (5.41) 

t:  z I 

and 


it  ° ''t  ®t 

where 


6^  is  defined  by 

dy^  2a  x^  + by^  + d 
tan  0^  - - bx^  + 2cy^  + e 


(5.42) 


(5.43) 


The  distance  from  the  missile  to  the  target  can  be  expressed 

in  terms  of  the  missile  position  coordinates  X„  and  ob- 

mm 

tained  by  solving  Eqs.  (5.29)  and  (5.30),  respectively, 
and  the  target  coordinates  x^  and  y^  as 
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(5.44) 


‘ <*t  - + ‘^t  - V"> 


The  miniimam  value  that  obtains  during  a particular  flight 
is  defined  as  the  "miss  distance”  (j>^,  l.e., 


4 = min(d^.^ 

m X r 


The  performance  of  the  missile  in  flight  depends  strongly 


(5.45) 


on  the  parameter  i|»^. 


5.6  Wind  Model 

-Ci 

To  simulate  the  effects  of  cross-winds  on  the  missile 

* 

trajectory,  a stochastic  wind  model  is  Incorporated  into 
CADLAW.  The  wind  is  assumed  to  blow  in  the  y-direction 
(see  Figure  5.1)  only  and  to  have  a magnitude  given  by 

Wy  = Wyjjj  ( 1 + 2o(R  - 0.5)]  (5.46) 

where  w is  the  mean  velocity  of  the  wind,  a is  the  standard 
yra 

deviation  of  the  distribution  about  the  mean,  and  R is  a 

random  number  between  zero  and  one . The  value  of  w ^ may  be 

ym  ■* 

varied  by  the  designer. 

53 


•’TTywo’-'^'jiyg ; 


6.  OPTIMIZATION  ALGORITHM 

Perhaps  the  key  element  of  the  CADLAW  system  is 
the  optimization  process  used  to  determine  the  "optimum 
design."  As  described  in  Section  2,  the  optimization 
process  involves  a pattern  search  of  a hypersurface,  defined 
by  an  objective  function,  to  obtain  a design  point  which  gives 
a maximum  or  minimum  value  for  the  function.  The  objective 
functions  employed  in  CADLAW  have  previously  been  defined  in 
Section  2.  The  obj«£ctive  of  this  section  is  to  describe 
the  particular  optimization  method  utilized  in  CADLAW. 

The  particular  alogrithm  adapted  to  the  computer 
program  CADLAW  is  the  flexible  tolerance  method  originally 
described  by  Paviani  and  Himmelblau  [12]  . The  remaining 
discussion  in  this  section  represent  mainly  excerpts  from 
Reference  13.  The  flexible  tolerance  method  can  be  classed 
as  a "search"  method  as  opposed  to  the  derivative- type  methods 
of  optimization.  In  the  purest  of  the  search  methods,  the 
directions  of  minimization  are  determined  solely  from  suc- 
cessive evaluation  of  the  objective  functions  f (x) . 

I As  a general  rule,  in  solving  unconstrained  nonlinear 
programming  problems,  gradient  and  second-derivative  methods 
converge  faster  than  direct  search  methods.  However,  in 
practice,  the  derivative-type  methods  have  two  main  barriers 
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to  their  implementation.  First,  in  problems  with  a, 
modestly  large  number  of  variables,  it  is  laborious  or 
impossible  to  provide  analytical  functions ' for  the  deriva- 
tives needed  in  a gradient  or  second-derivative  algorithm. 
Although  evaluation  of  the  derivatives  by  difference 
schemes  can  be  substituted  for  evaluation  of  the  analytical 
derivatives  the  numerical  error  introduced,  particularly 
in  the  vicinity  of  the  extremum  can  impair  the  use  of  such 
substitutions.  In  principle,  symbolic  manipulation  to 
evolve  analytical  derivatives  is  possible,  but  this 
technique  still  requires  considerable  development  bSfore 
it  becomes  a feasible  tool  in  practice.  ‘ In  any  case,  search 
methods  do  not  require  regularity  and  continuity  of  the 
objective  function  and  the  existence  of  derivatives. 

Second,  and  a related  point,  optimization  techniques  based 
on  the  evaluation  of  first  and  possibly  second  derivatives 
require  a relatively  large  amount  of  problem  preparation  by 
the  user  before  he  introduces  the  problem  into  the  algorithm, 
as  compared  with  search  techniques. 

Because  of  the  difficulties  described  above,  direct 
search  optimization  algorithms  have  been  devised  that, 
although  slower  to  execute  for  simple  problems,  in  practice 
may  prove  more  satisfactory  from  the  user's  viewpoint  than 
gradient  or  second-derivative  methods,  and  may  cost  less 

I 

to  use  if  the  cost  of  problem  preparation  time  is  high  rela- 
tive  to  the  computation  time. 
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The  general  nonlinear  progrcunining  problem  can  be 
stated  as  follows: 

Minimize:  f (x)  xeE^ 

Subject  to:  h. (x)  = 0 i = 1,  . . M (6-1) 

I 

g^(x)  ^Oi  = m+l,  . . .,p 

where  f (x) , h.  (x) , and  g . (x)  may  be  linear  and/or  nonlinear  ,, 

A#  X ^ X ^ 

functions.  In  many  nonlinear  programming  methods  a con- 
siderable portion  of  the  computation  time  is  spent  on 
satisfying  rather  rigorous  feasibility  requirements.  The 
flexible  tolerance  algorithm  [12],  on  thd  other  hand,  improves 
the  value  of  the  objective  function  by  using  information  pro- 
vided by  feasible  points,  as  well  as  certain  nonfeasible 
points  termed  near-feasible  points.  The  near-feasibility 

t 

limits  are  gradually  made  more  restrictive  as  the  search 

proceeds  toward  the  solution  of  the  programming  problem, 
until  in  the  limit  only  feasible  x vectors  in  Eos.  (6-1)  are 
accepted.  As  a result  of  this  basic  strategy  problem,  Eq.  (6-1) 

I , 

can  be  replaced  by  a simpler  problem,  having  the  same  solution: 
Minimize:  f (x)  xcE*^ 

Subject  to:  - T(x)  >0  (6-2) 

(k) 

where  4 ' is  the  value  of  the  flexible  tolerance  i qriterion 

for  feasibility  on  the  kth  stage  of  the  search  as  defined  by 
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Eq.  (6-3)  below,  and  T(x)  is  a positive  functional  of  all 
the  equality  and/or  inequality  constraints  of  Problem  (6-1) 
used  as  a measure  of  the  extent  of  contraiiit  violation. 

In  the  particular  application  here^  f(x)  is  the 
squared  difference  between  the  total  strain  energy  density 
and  the  element  strain  energy  density  summed  over  ^11  the 
elements  of  the  shell.  Hence,  f (x)  is  a nonlinear  function 
of  x.  The  h,  (x)  and  g . (x)  are  linear  constraint  relations 
for  the  composite  shell  optimization  problem. 


6.1  Flexible  Polyhedron  Search 

Nelder  and  Mead  [15  ] proposed  a method  of  search  which 
has  proved  to  be  an  effective  strategy  and  one  which  is 
easily  implemented  on  a digital  computer.  The  method  of 
Nelder  and  Mead  minimizes  a function  of  n independent 
variables  using  (n  -f  1)  vertices  of  a flexible  polyhedron 
in  E^.  Each  vertex  can  be  defined  by  a vector  x.  The 
vertex  (point)  in  E*^  which  yields  the  highest  value  of 
f (x)  is  projected  through  the  center. of  gravity  (centroid) 
of  the  remaining  vertices.  Improved  (lower)  values  of  the 
objective  function  are  found  by  successively  replacing  the 
point  with  the  highest  value  of  f (x)  by  better  points  until 
the  minimum  of  f (x)  is  found. 

The  details  of  the  algorithm  are  as  follows. 


JbVb  Xi  1^4  1 »•••»  *4  4 »»»*»*4.*  * t 


(k)  ^(k)  ^(k),T 

'***'*in 


rn  + 1, 


1 


• .(  ff.  'w-  u 


be  the  ith  ver1:ex  (point)  in  e”  on  the  )cth  stage  of  the 

search/  k » 0/  1/.../  and  let  the  value  of  the  objective 
fk^  Ik.) 

function  at  ' be  f(x.'  '),  In  addition#  We  need  to  label. 
X vectors  in  the  polyhedron  that  give  the  maximum  and 

mt 

minimum  values  of  f (x) . 

Define 

f(Xj^^^)  = max  {f  ) / . . . / 
with  the  corresponding  x|^^  = , and 


f ) ■ min  {f  )#.../  f J 


(k)  (k) 

with  the  corresponding  x.  * x ' , Since  the  polyhedron 

«wX 

in  e”  is  made  up  of  (n  +1)  vertices,  Xj^,,.,#  Xj^^_x» 

^n+2  centroid  of  all  the  vertices  excluding  Xj^. 

The  coordinates  of  the  centroid  are  given  by 


X 


(k) 

n+2,j 


1 

n 


n+1 
I(  I 
i-1 


j ■«  1/ . . . / n 


(6-3) 


where  the  index  j designates  each  coordinate  direction. 

The  initial  polyhedron  usually  is  selected  to  be  a 
regular  simplex  (it  does  not  have  to  be),  with  point  1 as 
the  origin,'  or  perhaps  the  centroid  as  the  origin,  as  in 
the  computer  code.  The  procedure  of  finding  a vertex  in 
e"  at  which  f(x)  has  a better  value  involves  four  operations: 
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i_-  w.-l.-.t* -B.C. - 


.,  Reflection.  Reflect  through  the  centroid  by 

computing 

where  a > 0 is  the  reflection  coefficient. 


x^^I  » centroid  computed  by  Eq.  (6“3) 
Cn+2 


- vertex  at  which  f(x)  is  the  largest  of 
<»h  " 

(n  + 1)  values  of  f(x)  on  kth  stage 


2.  Expaniion.  U < f (x<’‘> ),  expand  th.  v.etor 

<»n«  * 5n«>  •’y 


(k)  _ ^(k)  . - x^^h 

?n+4  *n+2  ^ ^^^n+3  ~n+2 


(6-5) 


where  y > 1 is  the  expansion  coefficient.  If 

replace.  by  x^+4  and  continue 
from  step  1 with  k - k + 1.  Otherwise,  replace  by 

and  continue  from  stop  1 with  k • k + 1. 

e»n+3 


Contractions 

. If 

*<?n+3>  ■ 

. fufh 

for  all  i 

h 

contract  the 

^ / (k) 

vector  (Xj^ 

computing 

x^^^  - 

^n+5 

^n+2 

- x*’'*) 

^n+2^ 

(6 

-6) 

where  0 < 3 < 1 is  the  contraction  coefficient.  Replace 


(k) 


*ti  ■ *>y  in+5 


and  return  to  step  1 to  continue  the 


search  on  the  (k  + l)st  stage. 


If 

*<!!n+3> 

> f(X^^^), 

,(k) 

n+3 

^h 

i • X f • * • t 

(k) 


from  X,  ' by  computing 


Si 


+ 0.5 (x!^^ 

••I  ; 


(6-7) 


and  return  to  step  1 to  continue  the  search  on  the 

»» 

(k  + l)st  stage. 

I 

The  criterion  used  by  Nelder  and  Head  to  terminate 
the  search  Was  to  test  to  determine  if 


, n+1 

T 

» + 1 i-1 


If(x{^^ 

A*  X 


- f(x^+2^l“  >'^<  e (6-8) 


(k) 

where  e la  an  arbitrarily  small  number,  and 

(kl 

value  of  the  objective  function  at  the  centroid  x^^2* 

6.2  The  Tolerance  Criterion 

The  tolerance  criterion  4 in  Eq.  6-2  is  selected  to 
be  a positive  decreasing  function  of  the  vertices  of  the 
flexible  polyhedron  in  e“; 

^r+1'  ^r+2^  * function  ♦ acts  as  a tolerance  criterion 

for  constraint  violation  throughout  the  entire  search, 
and  also  serves  as  a criterion  for  termination  of  the 
search.  Many  alternative  definitions  of  9 are  possible, 
but  the  one  Incorporated  into  the  algorithm  to  be  described 
is 


« 


(k) 


^ ^ i-1 


xf^> 

■•X 


„(k) 

Cr+2 


)♦ 


(0) 


- 2 (m  1}  t 


(6-9) 
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where 


t = size  of  initial  polyhedron 

m = number  of  equality  constraints 

' - ith  vertex  of  polyhedron  in  e" 

r **  (n  - m)  = number  of  degrees  of  freedom 

of  f (x)  in  Problem  (4-1) 

(k) 

X. i = vertex  corresponding  to  centroid  as  defined 
by  Eq.  (4-3),  with  n = r 
k » 0,  1,...  is  an  index  referring  to  number  of 
completed  stages  of  search 

(k-1) 

♦ = value  of  tolerance  criterion  on  (k  - l)st 

stage  of  search 
Let  the  second  term  in  the  braces  of  Eq.'  (6-9)  be 
denoted  by  0 . 


(k) 


m + 1 
r + 1 


r+1 

Z 

i=l 


-I 


x^^> 

Cr+2 


m + 1 
r + 1 


r+1 

Z 

i*l 


n 

Z 

j=l 


(Xij 


x^^^  .)2 

^r+2, 


(6-10) 


(k) 

where  » ii  “ 1,...,  n,  are  the  coordinates  of  the  ith 

vertex  of  the  flexible  polyhedron  in  E*^.  Observe  that  6^^^ 

(k) 

represents  the  average  distance  from  each  x.  , i « 1,..., 

1 

r + 1,  to  the  centroid  polyhedron  in  e”.  To 
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* (k) 

understand  the  behavior  of  4>  it  is  necessary  first  to 

understand  the  behavior  of  6.  It  is  obvious  that  the  value 

of  6 will  depend  on  the  size  of  the  polyhedron  in 

which  may  remain  unchanged,  expand,  or  contract,  depending 

on  which  one  of  the  four  operations  described  in  Sec.  4.1 

(k)  (k+1) 

is  used  to  carry  out  the  transition  from  x . to  Xj  ' . 

X *^x 

(k) 

Thus  4 'behaves  as  a positive  decreasing  function  of  x, 

*0 

(k) 

although  6 may  increase  or  decrease  during  the  progress 

of  the  search,  and  as  the  solution  of  the  problem  is  ap- 

(kl  (k) 

proached,  both  6 and  4 ' approach  zero 


(0) 


> ♦ 


(1) 


> 4 


(k) 


> 0 


(6-11) 


In  the  method  of  Nel^er  and  Mead,  when  it  is  not 
possible  to  find  better  values  of  f(x)  by  Eg.  (6-4^  the 

•m 

vertices  of  the  flexible  polyhedron  are  drawn  nearer  and 
nearer  to  that  vertex  corresponding  to  the  best  value  of 

I 

the  objective  function.  In  the  limit  complete  collapse  of 
all  the  vertices  of  the  flexible  polyhedron  takes  place 
onto  the  stationary  solution  of  f (x) . Thus,  as  the  search 
approaches  the'  stationary  solution  of  f (x) , the  value  of 

•m 

0^^^  given  by  Eg.' (6-10)  becomes  progressively  smaller 
because  the -average  distance  between  the  vertices  and  the 
centroid  of  the  polyhedron  shrinks  to  zero.  Since  on  each 

(k)  ' 

kth  stage  of  the  search  4'  'is  set  egual  to  the  smaller  value 
of  either  or  the  tolerance  criterion 


also  collapsed  and  in  the  limit. 


[ 


lim  = 0 

x-»-x*  (6-12) 

6. 3 Criterion  for  Constraint  Violation 

Consider  now  a functional  of  the  equality  and 
inequality  constraints  of  Problem  (4-1). 

m Pi. 

T(x)  ■ + I Z h.  (x)  + Z Ujg?(x)]*  (6-13) 

i=l  ^ " i=m+l  ^ ^ " 

i 

where  is  the  Heaviside  operator  such  that  ■ 0 ^or 

gi(x)  > 0 and  ■ 1 for  g^ (x)  < 0.  Therefore  T(x)  is 

defined  as  the  positive  square  root  of  the  sum  of  the 

squared  values  of  all  the  violated  equality  and/or  inequality 

constraints  of  Problem  (6-1).  Note  that  T(x)  > 0 for  all 
n 

xeE  . In  particular,  if  Z h.  (x)  is  convex  and  the  g^  (x) , 

i=l  1 

i * m 4-  1, . . . , p,  are  concave  functions,  then  T(x)  is  a convex 
function  with  a global  minimum  T(x)  0 for  all  feasible 
X vectors;  i.e.,  for  any  {xjh. (x)  = 0,g. (x)>0  for  i * l,...,p}. 
Also,  T(x)  >0  for  all  x vectors  that  are  nonfeasible.  For 
a given  x^^^es”,  the  value  of  tCx)  evaluated  at  x^^^  using 
Eq.  (4-13)  can  be  used  to  distinguish  between  feasible  and 
nonfeasible. points.  If  T(x^^^)  = 0,  x^^^  is  feasible; 
if  T(x^^^)  > 0,  x^^^  is  nonfeasible.  On  the.  other  hand,  a 
small  value  of  T(x^^^)  implies  that  x^^^  is  relatively  near' 

00  0t 
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' (J,) 

to  the  feasible  region,  and  a large  value  for  T(x  ) 

Oc) 

implies  that  x is  relatively  far  from  the  feasible  region. 

6 . 4 Concept  of  Near-Feasibility 

Near-feasible  x vectors  are  those  points  in 
•>0 

that  are  not  feasible,  but  nevertheless  =ilmost  feasible, 
in  the  sense  given  below.  To  establish  clear-cut  distinc- 
tion between  feasible,  near-feasible,  emd  nonfeasible  points, 
let  4 be  the  value  of  4 on  the  kth  stage  of  the  optimization 
search  and  let  x^^^  be  any  vector  in  The  vector 

• ' a# 

is  said  to  be 

1.  Feasible,  if  T(x^^^)  * 0 

2.  Near-feasible,  if  0 < T(x^^^)< 

3.  Nonfeasible,  if  T(x^^^)  > 

Thus  the  region  of  near-feasibility  is  defined  as 

♦ - T(x)  > 0 (6-14) 

(k)  (k+1) 

On  any  transition  from  x to  x , the  move  is  said  to 
be  feasible  if  T(x^*^^^^)  « 0‘,  near-feasible  if 

0 < T(x^^‘*’^M  < and  nonfeasible  if  T(x  . 

Note  that  the  value  of  0 on  the  (k  + l)th  stage  of  the 
search  is  determined  only  after  has  been  located  as 

I 

either  a feasible  or  ‘feasible  point. 
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6,5  Strategy  of  the  Flexible  Tolerance  Algorithm 

In  this  section  it  is  demonstrated  that  the  general 
nonlinear  programming  Problem  (4-1)  can  be  replaced  by  the 
easier  problem  of  minimizing  f(x)  subject  to  one  gross 
inequality  constraint  as  follows: 


Minimize: 

f (X) 

xeE^ 

Subject  to: 

- T(x)  > 0 

(6-15) 

The  flexible  polyhedron  search  of  Nelder  and  Mead  is  a 
convenient  and  effeo.'yive  but  not  essential  method  of  minimiz- 
ing f (x)  as  an  unconstrained  function  when  the  constraint 
in  (6-15)  is  not  active,  and  is  also  used  to  minimize  T(x)| 

to  satisfy  the  single  constraint  in  (6-15)  when  the  con- 

(k) 

straint  is  active.  The  general  strategy  is  to  reduce  v 
as  the  search  progresses,  thus  tightening  the  region  of  near- 
feasibility, and  to  segregate  the  minimization  of  f (x)  from 
the  steps  taken  to  satisfy  the  constraint  in  (6-15)  . For 
a given  value  of  the  value  for.T(x)  at  will  be 

either  (1)  T(x^^"*^^^)  < in  which  case  is  either 

a feasible  or  a near-feasible  point  and  will  be  accepted  as 

a permitted  move,  or  (2)  T(x^*''*’^^)  > 4^^^,  in  which  case 

^(k+1)  classed  as  nonfeasible,  and  an  x vector  closer  to 

or  in  the  feasible  region  must  be  found  in  lieu  of 

A* 

One  way  of  getting  an  closer  to  the  feasible  iegion  i^ 

to  minimize  the  value  of  T(x^^^^^)  as  defined  by  Eq.  (6-13) 
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until 

To  demonstrate  that  the  solution  of  Problem  (6-15 

is  equivalent  to  the  solution  of  Problem  (6-1) , it  is 

(k) 

sufficient  to  consider  the  behavior  of  ^ . Because 

(kl 

4 ' is  a positive  nonincreasing  function  such  that 

a 0 only  when  it  is  no  longer  possible  to  Improve  the 

value  of  f (x)  in  Problem  (6-15)  the  region  of  near-feasibility,' 

given  by  Eq.  (6-14)  is  gradually  restricted  as  the  ^earch 

proceeds  toward  the  solution  of  Problem  (6-15) . In  the 

(k) 

limit/  that  is,  when  all  the  vertices,  x.  , i ■■  1,..., 

r.-f  1/  of  the  flexible  polyhedron  in  have  collapsed  into 

one  single  point  at  x*,  then  « 0 and  only  x vectors  that 

are  feasible,  that  is,  (x  jhj^(x)  » 0,  9j^(x)  > 0 for  i}, 

can  satisfy  the  requirements  of  the  inequality  in  Eq.  (4-14). 

(k) 

In  other  words,  if  4'  ' » 0,  since  T(x)  cannot  be  negative, 
the  only  possible  value  for  T(x)  is  T(x)  ■ 0,  which  requires 

I 

that  all  the  constraints  of  Problei  (6-1)  . be  satisfied. 

t 

Because  the  tolerance  criterion  4 is  a positive 
nonincreasing  function  of  the  sequence  of  points 

m0 

x^^^,...,  X*  generated  during  the  progression  of  the  search, 
because  4 does -not  depend  on  the '.value  of  the  objective 
function  nor  on  the  values  of  the  constraints,  and  because 
in  the  limit,  4*  ■ 0,  convergence  of  the  algorithm  is  as- 
sured for  the  following  reasons  i 


1.  The  manner  in  which  ^ is  computed  by  Eq.  (6-9) 
prevents  the  tolerance  criterion  from  increasing. 

J 

If  $ were  allowed  to  increase  without  bound,  the 
possibility  would  arise  of  being  able  to  improve 
the  value  of  f(x)  at  the  expense  of  getting 
further  and  further  away  from  the  feasible  region. 

2.  When  the  optimal  solution  of  Problem  (6-1)  is  an 
interior  point  (no  equality  constraints) , con- 
vergence of  the  algorithm  is  assured  because  of 
the  property  of  the  flexible  polyhedron  of  col- 
lapsing  only  when  approaching  the  optimum  of 

of  f(x)  in  Problem  (6-1).  in  these  circum- 
stances  T(x)  has  no  effect  on  the  convergence 
of  the  algorithm  because,  in  the  final  stages  of 

If 

the  search,  x^,  i * 1,...,  r + 1,  are  interior 

/If  \ 

points  yielding  a T(x'  ')  » 0,  which  implies 

•V 

that  inequality  (6-14)  is  satisfied  for  all 
(k) 

X.  and  that  Problem  (6-1)  has  not  active  con- 
straints. 

3.  When  the  optimum  of  Problem  (6-1)  is  not  an 
interior  point  [either  because  x*  is  a boundary 
point  or  because  Problem  (6-i)  includes  only 

' equality  constraints] , convergence  is  assured 

t <• 

because  of  the  condition  imposed  by  inequality 
(6-14 )r  that  is  - T(x)  > 0.  The  flexible 
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poLS'hedron  will  not  collapse  as  long  as  it  is 

possible  to  find  a better  f ) such  that 

♦ > 0. 

(k) 

Let  'be  the  vertex  of  the  flexible  polyhedron 

such  that  - T(xf^^)  > 0;  f(xi^^)  is  the  best  viluo  o£ 

f(x)  obtained  on  the  kth  stage  of  the  search.  Let 
~ ««x 

fk) 

be  the  vertices  obtained  during  any  reflection  of  ' 

through  the  centroid  of  the  polyhedron  such  that 

- T(x|^‘*‘^^)  >0.  If  > f(Xj[^^  for  every 

(k) 

reflection  of  x^  through  the  centroid,  the  values  of 

will  decrease  because  of  contractions  in  the  flexible 

/k) 

polyhedron.  In  such  a case,  the  values  of  'are  reduced, 

(k) 

and  Xj^  roust  satisfy  the  constraints  of  Problem  (4^1)  nore 
and  more  closely  until  the  search  is  terminated  because 
♦ ^ e. 

On  the  other  hand,  if  f(x5^’*’^^)  < f(xf^^),  the  value 

. ” <•*> 

(k) 

of  4 ' ' will  not  decrease  because  no  contraction  of  the 
polyhedron  takes  place  and  is  replaced  by  a better 

vertex.  As  long  as  it  is  possible  to  determine  either  a 
feasible  or  near-feasible  point  such  that  f ) < f } , 

there  will  be  reflections  and  expansions  of  the  flexible 
polyhedron.  Thus  premature  terynination  of  the  search  at  a 
nonlocal  optimum  is  avoided  because  the  polyhedron  will  not. 
collapse  if  there  exists  an  such  that 

"*X 

^ T(xf^'*‘^^)  > 0 and  fixi^'*’^^)  < f(x*  4>  s).  (6-15) 

^ * •“  ••  as  • 

/ 
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One  advantage  of  the  flexible  tolerance  strategy 

1 

typified  by  Problem  (6-15)  is  that  the  extent  of  the 
violation  of  the  constraints  included  in  Problem  (6-1)  is 
progressively  decreased  as  the  search  moves  toward  the 
solution  of  Problem  (4-1).  Because  the  equality  and/or  the 
inequality  constraints  are  loosely  satisfied  in  the  early 
stages  of  the  search,  and  more  tightly  satisfied  only  as  the 
search  approaches  the  solution  of  Problem  (6-15),  the 
overall  compution  effort  required  in  the  optimization  is 
cohsiderably  reduced. 

Another  advantage  of  the  flexible  tolerance  strategy 

I 

Oc) 

is  that  ♦ can  be  conveniently  used  as  a criterion  for 

termination  of  the  search.  For  all  practical  purposes  it  is 

(k) 

sufficient  to  continue  the  search  until  $ becomes  smaller 

i 

than  some  arbitrarily  selected  positive  number  c.  In  the 

(k) 

final  stages  of  the  search,  4>  is  also  a measure  of  the 

(k) 

average  distance  from  each  vertex  , i » 1,...,  r + 1, 

to  the  centroid  polyhedron  in  e”.  If  < e, 

(k) 

a substantial  number  of  the  vertices  x.  are  contained  in 

X 

a hypersphere  of  radius  e.  (If  the  last  polyhedron  of  the 

(k)  ' 

search  were  regular,  all  the  vertices  x.  would  be  contained 

X 

by  the  hyper  sphere  of  radius  c,  but  because  the  polyhedron 

is  distorted,  some  vertices  may  be  outside  the  sphere.) 

(k)  ' 

Therefore,  if  « ' < c,  the  chances  are  that  the  value  of 

■ I ■ I 

f (x)  cannot  be  improved  without  having  t6  further  reduce  the 

m0 
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size  of  the  polyhedron.  This  implies  that  a change  of 
(k) 

2e  in  the  x.  corresponding  to  the  best  value  of  f (x) , 

— X ~ 

(k) 

that  is,  x^  , will  not  improve  the  value  of  the  objective 
function.  Hence,  upon  termination  of  the  search,  the 
following  condition  is  satisfied: 


f(Xj[^^  < f(x*  + e ) 


(6-16) 


(k) 

Since  Eq.  (4-14)  is  satisfied  for  every  move,  if  ♦ £ e, 

(k) 

it  IS  obvious  that  the  condition  e - T(x^  ) ^0  is  also 


satisfied,  or 


(k) 


m 


1 


T(Xj^  ')  ■ [ Z h?(x)  + Z U . (x)g.  (x)  ]^  < e (f-17) 

i=l  ■“  i=m+l 


Equation  (6-17)  implies  that,  upon  termination  of  the 
search,  the  combined  value  of  all  the  violated  constraints 
does  not  exceed  e.  Certainly,  no  individual  constraint  can 
be  violated  by  more  than  c either. 

Thus , by  conducting  the  search  in  the  manner  described 
one  arrives  at  an  "optimum  design"  point.  This  procedure  is 
employed  in  CADLAW. 
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7.  CADLAW  DESCRIPTION  AND  USER  INSTRUCTIONS 


The  design  program  CADLAW  operates  interactively  from 
a demand  terminal.  The  program  has  the  capability  to  perform 
optimized  design  studies  or  classical  parameter  studies  and 
to  present  the  results  in  either  a tabular  or  graphical  form. 

To  operate  the  program  a basic  awareness  of  the  organization 
of  the  program  is  helpful. 

7. i Design  Program  Organization 

As  described  earlier  in  Section  2,  CADLAW  is  designed 

to  permit  study  of  three  different,  independent  subjects  re- 

I 

lated  to  light  antitank  weapon  systems.  The  first  subject 
concerns  the  propulsion  system,  including  structural  features,  » 

of  a conventional  in-tube  burning  rocket  system.  The  second 
subject  concerns  the  propulsion  system,  including  structural 
features,  of  a recoilless,  closed  rocket  system.  The  third 
subject  with  which  CADLAW  deals  is  the  flight  or  trajectory 
phase  of  operation  of  a light  antitank  weapon  system.  The 
different  type  propulsion  systems  can  obviously  be  studied 

T 

independently  since  most  antitank  systems  will  have  only  one 
propulsion  system.  The  trajectory  phase  of  antitank  weapon 

i 

operation  can  be  studied  independent  of  the  propulsion  phase 
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since  all  propulsion  is  accomplished  within  the  launch  tube 
before  flight  characteristics  become  of  interest.  The  CADLAW 
user  must  decide  which  of  the  three  subject  areas  is  to  be 
studied  and  direct  the  computer  to  implement  the  subject 
matter  of  interest.  Each  of  these  three  subjects  are  contained 
in  three  separate  overlays  and  are  loaded  and  unloaded  from 
the  computer  core  as  the  user  directs.  The  basic  overlay 
design  of  CADLAW  is  described  schematically  in  Table  7.1. 

Within  each  of  the  basic  overlays  dealing  with  one  of 
the  three  subjects  described  above,  there  exists  a capability 
to  either  perform  a^  optimized  design  study  or  classical  para- 
meter study  as  the  user  desires.  Only  one  type  of  Study  may 

; 

be  performed  lU  a time.  However,  both  type  studies ‘may  be 
performed  by  operating  in  a sequential  mode.  These  components 
of  the  program  are  also  illustrated  in  Table  7.1. 


7.2  User  Instructions 

The  CADLAW  system  is  basically  semi-tutorial  in  nature, 
i.e.,  if  desired,  the  system  explains  the  basic  organization 

I 

of  itself  and  informs  the  user  of  the  responses  to  be  made  to 
implement  the  system.  To  save  time  the  system  will  skip  over 
all  instructions  if  the  user  is  ’familiar  with  it. 

Since  the  CADLAW  system  is  interactive,  the  system  com- 
municates with  the  user  by  printing  messages  on  the  CRT 
screen  of  the  interactive  terminal.  These  messages  are  either 
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Table  7.1  Description  of  CADIAW  Overlay  Structure 


instructions  on  program  usage,  questions  for  the  user,  or 
statements  of  results  or  data.  The  user  responds  to  the 
computer  by  typing  information  on  the  terminal  keyboard. 

Since  the  program  is  self  explanatory,  no  specific  instruc- 
tions need  be  supplied  the  user.  Only  those  data  which  differ 
from  a basic  data  set  defined  within  the  system  required  in- 
put. All  others  assume  default  values. 

7.3  Sample  CADLAW  Execution 

It  appears  that  familiarity  with  the  program  car,  best 
be’ gained  by  following  through  a typical  execution  cycle.  The 
following  pages  represent  copies  of  the  CRT  screen  displays 
as  seen  by  the  user.  When  a right  bracket  appears  in  column 
one,  the  computer  is  asking  for  a response  from  the  user. 

All  data  input  to  the  program  is  in  a free-format  style.  Hence 
the  user  may  simply  type  in  the  desired  pieces  of  information, 
each  separated  by  a comma.  A read  sequence  is  terminated  by 
typing  in  a series  of  zeros  for  the  information  asked  for. 
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